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Abstract 

We solve a 2D model of A^-component dense electron gas in the limit N ^ oo and in a range of 
the Coulomb interaction parameter: N^^^"^ ^ ^ 1. The quasiparticle interaction on the Fermi 
circle vanishes as /Nm. The ground state energy and the effective mass are found as series in 

2 /3 

powers of . In the quantum Hall state on the lowest Landau level at integer filling: 1 ^ < A^, 
the charge activation energy gap and the exchange constant are: A = \og{rsN^/'^)fiuj}j /i^ and 
J = Q.QQhuJH /t-'- 

PACS numbers: 71.10.Ca, 73.43. Cd 
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Two-dimensional electron gas (2DEG) in GaAs quantum wells and Si heterostructures 
|1|] is a unique system where density of electrons n can be varied widely. Effects of 
electrostatic Coulomb interaction are determined by the variable dimensionless coupling: 
Tg = e^m/y/Trnfi^. In Si MOSFET interesting phenomena was observed recently at rela- 
tively large rj ~ 10: a conductivity 'fan' similar to the metal-insulator transition 0, a 
sharp increase of the effective electron mass and magnetic susceptibility At large no 
exact model exists and many phenomenological approximations have been developed. 

At small Tg there is an exact model - dense electron gas 0] but two of its key predictions 
have not been confirmed in the experiment. First, theory predicts a reduction of the effective 
mass: m*/m— 1 ~ log(l/rs) at ^ - because the forward scattering is larger then the 
backward scattering for Coulomb potential. Second, in the integer quantum Hall ferromagnet 
state on the lowest Landau level, theory predicts the Coulomb charge activation gap A ~ 
e^^eH/hc in strong magnetic field H (r, -> 0). Yet, in Si(lOO) MOSFET the Shubnikov-de 
Haas oscillation experiments 0, find: m* /m = 1 + COSr^, at 3 > > 0.9, and the 
magnetocapacitance experiment [§] finds much lower charge activation gaps: A ~ heH/mc, 
at 4 > > 1.5. In phenomenological models (see Refs. in [|^) experimental m*{rs) > m 
was reproduced in the large domain with a crossover at r* ~ 1 from the theory with 
m* < m to the large with m* > m. Experiments bound this crossover to be small 
r* < 0.9. But perturbation series in one parameter could hardly have the two so widely 
different crossovers: r] and r*. 

In this Letter we show that a systematic model of multicomponent dense electron gas 
agrees with the experiment. In Si electron states have a valley degeneracy |l|, corresponding 
to different band minima. For (100) orientation of 2D plane there is iV = 4 equivalent 
spin- valley states. Two valley band states here differ only by plane- waves exp(±zQz) in the 
perpendicular direction with an atomic value of momentum Q. These states are strongly 
orthogonal. For (111) orientation the spin-valley degeneracy is especially large = 12, 
but in experiments only A^ = 4 states are filled with electrons ll|] - a long standing puzzle. 
Electron gas in the limit N ^ oo was first considered in Ref . [^] and the leading term of the 
ground state energy was shown to coincide with that of a dense charged boson gas |T0[. 

For A^-component dense 2DEG in the limit A^ — oo we find a crossover at rl = l/A^3/2 
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from the theory at vanishing to a novel systematic theory at: 

< 1 and r, > r* = l/N^'"^. (1) 

It describes a system of interacting electrons and plasmons. Plasmons are excitations with 
large characteristic momentum qo ^ Pf and energy ujq ^ ep. Landau theory of Fermi liquid 
uses only one function - the amplitude of quasiparticles scattering - to derive all Landau 
Fermi liquid parameters. In our theory the interaction between quasiparticles vanishes as 
1/Nm and the electron subsystem is the ideal Fermi gas. The exchange of high energy 
plasmons gives rise to a polaronic like renormalization of the effective electron mass as a 
series in powers of r^/^. In the case of integer quantum Hall state on the lowest Landau level 
we predict a linear dependence on magnetic field of the charge activation energy gap and 
the exchange constant that agree better with the experiments. 



2DEG Hamiltonian is expressed in terms of the second quantized electron operators [|1T 
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H = — 
2m 

+ lJJ j^^PUx)4{r)M^MS) d'x d'r, (2) 

where a,f3 = 1..N are spin- valley indices. is an even length of fermion spinor. No valley 
mixing in the density operators, isotropic mass tensor and a positive uniform charge at 
large distance d from the 2D plane are assumed. At first we assume zero magnetic field 
and ap = e'^/hvp ~ TsVN -C 1 instead of the first condition (|I]) and then we extend the 
theory to all -C 1. We use Matsubara diagrammatic expansion in terms of the Coulomb 
interaction. New Coulomb line multiplies the diagram by a small parameter ap, whereas 
each electron loops brings a large factor with it. Therefore, the leading order is given by 
diagrams with minimum number of Coulomb lines per electron loop, with an essential block 
being the RPA screened Coulomb interaction - a sequence of alternating Coulomb lines and 
electron loops |]T2[. The second condition (|I]) makes a typical RPA momentum to be large: 



g ~ go ^ Pf, where go = (Svre^mn)^/^ = r^^y^pp and luq = c^l\f2m defines a typical 
plasmon momentum and energy. The electron polarization operator is: 

nw,g) = ^— , (3) 

where e(g) = q^/2m is the electron dispersion. It depends on the total density n and is 
independent of an electron distribution over A^ spin- valleys as long as pp^ <^ go- The RPA 




FIG. 1: a) Four-leg plasmon vertex and b) two three- leg plasmon vertices contribution to the 
energy. Solid lines are electron propagations and wavy lines are plasmon propagations. 



Coulomb interaction is: 



2716' 



-1 



(4) 



Its pole corresponds to a plasmon excitation with a dispersion: u{q) = \J q'^ + qlql2m. 
Electrons can be integrated out from the model (0) leaving an effective theory of plasmons 
with typical momentum go and a propagator (^) which weakly interact (in the limit ^ 1) 
due to three plasmon, four plasmon etc vertices. A plasmon vertex with k leg with external 
momenta qi ^ q and frequencies \uji\ ~ vanishes in the limit q]^ /m <^ \uji\ as nq^ /muj^ for 
k even and nq^ jrr^LO^^^ for k odd due to the electron number conservation. 

The kinetic energy density E^^m = im'^/Nm is small because electrons are evenly dis- 
tributed over spin-valley states. In the leading order the energy of 2DEG is given by 
the zero energy of plasmons: Eq = Ehf + E' , where Ehf = —Srsu'^/Smy/N is the 2D 
Hartree-Fock exchange energy and E' is the RPA energy [l^ : 

du d^q 



E'-- 



, , , 27re2 , ^\ 27re2 , _^ 
log ( 1 + q) n(^, q) 
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At g ~ Pi? we single out from E' a term that cancels E^p exactly and the remaining plasmon 
energy (also called a correlation energy) can be evaluated as [ll3| : 



= - 



m 



(6) 



Plasmon energy (|^) coincides in the leading order with the ground state energy of a dense 
charged boson gas [|1^]. In the case of anisotropic electron mass tensor of Si(lll) we compare 
two states with i) all A^ = 12 states being filled and ii) only A^ = 4 states related by the spin 
and time reversal being filled. Because log in (|^) is a convex function the plasmon energy 
(^) in the case ii) is always lower then in the case i) by an amount ~ 0.04r^/^n^/ (m). This 
can overcome the kinetic energy difference: 7rn^/6(m) at experimental densities > 1. 
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FIG. 2: A particle-hole ladder for vertex T. Two leg element D three leg element four leg 
elements r4 is shown. 

The next corrections to the plasmon energy are given by the second order diagram a) 
and the third order diagram b) on Fig.[l|. Provided the internal frequency and momentum 
are related by plasmon dispersion, the value of a diagram is proportional to r^/'^ in power of 
Ni — NL + 2, where Ni is the number of Coulomb lines and is the number of electron loops. 
We have evaluated numerically two diagrams on Fig.|T| and find: Eq + Ei = —(2.03191 r^^^ — 
0.156(l)r^) n'^/'m, where the first term represents Eq.(|]). Comparing vertices of the fermion 
(I) and the dense charged bose gas models we conclude that these models are different. 

A specific feature of our mean field theory that distinguishes it from the standard Landau 
Fermi liquid theory is that electron and plasmon excitations are found in a wide momenta 



range pp ^ p ^ qo- Let us prove that the four fermion vertex ||TT| vanishes on the infrared 
side of this range: Taa-i3i3{^i,Pi) ~ Q'^/mn at Pi/qo 0, where ~ pf/2m and i = 1..4. 
Q = ei — €2 and Q = Pi—p2 are the transfered frequency and momentum in the particle- hole 
channel on Fig.|^. A pair of electron and hole Green functions with independent integration 
over frequency and momentum and sum over spin- valley index 7 gives a large term |p!T| . 
Therefore vertex T{Q, Q) is a ladder of alternating pairs - polarization operators Q) - 
and blocks M, defined as a set of all diagrams that can not be cut over the pair lines. The 
block M is assumed to be independent of the internal integration momenta and frequency 
of the adjacent n(fi, Q). The block M can be further divided into i) a part that can be cut 
over one plasmon line: T^D{VL, (5)r3, where is a three leg one-plasmon irreducible part at 
vanishing leg momenta, and ii) a four leg one-plasmon irreducible part r4 at vanishing leg 
momenta. Let r£)(i7, Q) and rr(^^, Q) be the parts of the total vertex r(fi, Q) that end on 
the left with a plasmon line D or with a block The Dyson equation for the particle-hole 
channel is algebraic: 

/ 





TiD{Vl - 1) DV^Ii 

r4nr3 r4n 
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(7) 
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FIG. 3: Four-leg vertex is a sum of two diagrams in the leading order. 



The total four fermion vertex is given as F = F3F/) + Fp: 

^ DTl + T,{l-DU) 

Fi + (i + F4n)(i-Dn)- ^ ^ 

If ~ Q'^/m then we estimate the plasmon propagator D{Q,Q) ^ Il^^{Q,Q) ~ to be 
small but the factor: (1 — DU) ~ Q^, is even smaller in the limit Q ^ 0. We use the 
Ward identity: F3 = Z^^, where Z < 1 is the Green function pole renormalization. F4 is 
given in the lowest order by the two diagrams on Fig.i Tl = -v/¥F(2/3)F(5/6)rf V^Sm. 
It represents a retarded attractive interaction because it appears in the second order of 
perturbation in plasmon exchange. Assuming that F4 is finite we conclude from (§) that 
F(0,Q) = D{0,Q) = Qy4m*nZ\ 

This suggests an analogy between Fermi gas and a ferromagnet with spontaneously broken 
symmetry where both spin-wave dispersion and interaction vanish as according to the 
Goldstone theorem. A continuation to momenta on the Fermi circle gives the quasiparticle 
interaction: F = 2T[/mN . 

There is no infrared divergences in the block F4 despite of the possibility that some inter- 
nal frequency integration are determined by electron energy poles rather then by plasmon 
energy poles in the momentum range: pp '€i p <^ go- Would be divergent are those diagrams 
that i) can be cut into two or more parts over only electron lines; ii) each part, emerged as a 
result, is one-plasmon irreducible after all internal electron lines are contracted into points. 
These parts where all plasmon lines are closed into loops are fermion vertices symbolically 
written as J ip'^^ip^dx^dt. In the Cooper channel the four-leg vertex V4 = F ~ van- 



ishes thus cutting off a specific 2D logarithmic divergence |Tl[. All higher order vertices are 
regular as functions of leg momenta and are irrelevant in 2D as a simple power counting 
shows: t ~ a;^ and ip ~ 1/x. 

Despite the vanishing interaction between electrons there is an essential plasmon polaronic 
effect. The electron self-energy ||1 1|| : 



S(6, p) = -j Dicu, q)G{e + ^, p + q) (9) 



FIG. 4: Two diagrams for a spin susceptibility correction 



is related to the electron Green function: G~^(e, p) = ie — /2m + fi — S(e, p), where fi is 
the chemical potential. At e ^ cjo and p <^ go, we calculate the self-energy (H) and find: 
G^^i^jP) = Z {ie — e{p))^^ , where the quasiparticle dispersion is: e(p) = {p"^ — pp) / 2m* , with 
the renormalized effective mass: 



A static screened potential of an external charge z immersed into 2DEG is given by 
D{0,q) (^: V{q) = 2Tie^zq^ / {q^ + q^). It grows with the transfered momentum that helps 
to explain the heavier effective mass ([T0|) because the backward scattering is larger than the 
forward scattering. 

The temperature dependence of the effective mass comes from the momentum range: 
Pf ^ V ^ Taking into account the temperature dependence of the polarization operator 
in the equation for the self-energy we encounter logarithmic infrared divergence. Deriving 
and solving a simple renormalization group equation we find a growing effective mass if 
momentum of the quasiparticle decreases. On the Fermi line we find: 



Si(lOO) MOSFET the effective mass (|TT|) could strongly depend on temperature. 

At T <^ cuo plasmon thermal fluctuations freeze out whereas the ideal electron subsystem 
evolves from a degenerate gas at T ^ ei? to the Boltzmann gas at T ^ tp- 2DEG ther- 
modynamic potential: ^{T) = NQo{T,fi) + Eq + Ei, where Qo{T,fi) is the thermodynamic 
potential of the ideal Fermi gas with a property: NQq{0, ep) = E^in- 

Static magnetic susceptibility is related to a change of the total energy in external 
magnetic field: 6E = —Nx*{iJ'bHY/2. x* is given by the Pauli ideal gas susceptibility: 
X = m/27r, and by two diagrams on Fig.^. Evaluating them we find: x* Ix = m* /m. There 
is no exchange part coming from the interaction between quasiparticles. Thus, there is a 




(10) 



and Z ^ 



1 + r(l/3)r(7/6) r1^^ /2^/tx is a renormalization of the Green function pole. 
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single Landau Fermi liquid parameter m* /m given by a series in powers of r^/^. We ob- 
serve that the experimental effective mass |^ agrees well with (|10D- The experimental 
susceptibility 0, Q is larger then x* fhat indicates some exchange effects for = 4. 



In the end we consider the quantum Hall state |T5|, with an integer filling factor i/, 
defined as a ratio of electron number to the magnetic flux number. \i 1 <^ v < N then v is 
a number of occupied spin-valley states on the lowest Landau level. This state is degenerate 
under global rotations in electron Fock space by unitary matrices from the Grassmanian 
coset: F = U{N)/U{v) ® U{N — u), that corresponds to a ferromagnetic order parameter 
at zero temperature. We use the so-called magnetic units: h = 1, e = c, magnetic length 
Ih = cyclotron frequency ujh = eH/mc = 1/m and = y^e^ /coH^Hy/i^- The polarization 
operator is: 

and at qlu ^ 1 it transforms into @. We find the self-energy using the Green function: 
G^^(e) = ie + fi, on the lowest Landau level. The difference between the self-energies 
for electron in the empty spin-valley state and electron in the occupied spin-valley state 
gives the charge activation gap: A = Eg — Eq. We single out the Hartree-Fock term: 



Ehf = Y7r/2e and find: 

where only the real part of electron Green function is essential, fi = A/2 is the chemical 
potential inside the ferromagnetic gap and the last fraction in (|TB|) becomes approximately 
delta-function 6{uj), if /i <^ ujh- Hartree-Fock term cancels out exactly like in zero magnetic 
field and after frequency integration we find: 

A = ^ r D{q)dq = ^ (log(r.z.^/^) + 0.27?) , (14) 
Z Jo v ^ ' 

where L'-^(g) = e^'/Vr.v^ + z/E^i g^"~V2's!s. Note that A < Euf- The gap on the 

lowest Landau level ([l^) is similar to the gap in the N = 2 quasi-classical case of weak 



magnetic field and odd integer ^ 1 ||16 



We find a dispersion of spin waves as a pole of correlation function C{x — r) = 
{iljl{x)tlJo{x)'ipl{r)ipe{r)) , which is uniform because it describes a neutral excitation. It is 
given by a ladder set of diagrams like the second one on Fig.^, where electron Green func- 
tion include the self-energy: = ie ± A/2. Using Landau gauge representation of the 
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density operators we find: 

00/ j2 ^\ ^ 

a sum over A;- leg ladder diagrams. It has a pole at iQ = t{Q) with the spin- wave dispersion: 

<Q) = ^ r D{q) (1 - Jo(gg)) dq. (16) 

At intermediate wavelengths rgV^I'^ <^ QIh ^ 1 spin wave dispersion is logarithmic: t{Q) = 
—ujh ^og^Qlff) I whereas in the long wavelength limit at Ql^ ^ r^z/^/^ we recover Goldstone 
dispersion: e(Q) = JQ^, where the exchange constant can be evaluated numerically: 

J='^ r q^D{q)dq = 0.6613 ^ (17) 
Jo V 



A pair of Skyrmion topological defects of ferromagnetic order has a lower activation 
energy: A = J , then the electron-hole pair (0). 
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